Accepted by ApJL 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



X-RAY SPECTRAL VARIATIONS IN THE YOUNGEST GALACTIC SUPERNOVA REMNANT Gl. 9+0.3 

Stephen P. Reynolds, 1 Kazimierz J. Borkowski, 1 David A. Green, 2 Una Hwang, 3 Ilana Harrus, 3 & Robert 

Petre 3 

Accepted by ApJL 

ABSTRACT 

The discovery of the youngest Galactic supernova remnant (SNR) Gl.9+0.3 has allowed a look at 
a stage of SNR evolution never before observed. We analyze the 50 ks Chandra observation with 
particular regard to spectral variations. The very high column density (Njj ~ 6 x 10 22 cm -2 ) implies 
that dust scattering is important, and we use a simple scattering model in our spectral analysis. 
The integrated X-ray spectrum of Gl.9+0.3 is well described by synchrotron emission from a power- 
law electron distribution with an exponential cutoff. Using our measured radio flux and including 
scattering effects, we find a rolloff frequency of 5.4(3.0, 10.2) x 10 17 Hz (hv ro \\ — 2.2 keV). Including 
scattering in a two-region model gives lower values of u m \\ by over a factor of 2. Dividing Gl.9+0.3 
into six regions, we find a systematic pattern in which spectra are hardest (highest ^ ro ii) in the bright 
SE and NW limbs of the shell. They steepen as one moves around the shell or into the interior. 
The extensions beyond the bright parts of the shell have the hardest spectra of all. We interpret the 
results in terms of dependence of shock acceleration properties on the obliquity angle #Bn between the 
shock velocity and a fairly uniform upstream magnetic field. This interpretation probably requires a 
Type la event. If electron acceleration is limited by synchrotron losses, the spectral variations require 
obliquity-dependence of the acceleration rate independent of the magnetic-field strength. 

Subject headings: supernova remnants — ISM: individual (Gl.9+0.3) — X-rays: ISM 




1. INTRODUCTION 



The supernova remnant (SNR) Gl.9+0.3 has recently 
been shown to have expanded by about 16% between 
1985 and 2007, implying an age of order 10 years - the 



younges t supernova remnant in the Galaxy ( Reynolds et 
al. 2008 hereafter Paper I) . The expansion was confirmed 
with new VLA observations from March 2008 (Green et 
al.|2008^. The X-ray and current radio images are shown 
in r'ig.Jl] The X-ray spectrum is featureless and well de- 
scribeclby the loss-steepened tail of the synchrotron spec- 
trum inferred from radio frequencies. No thermal X-ray 
emission is apparent. Gl.9+0.3 provides a unique oppor- 
tunity to study a SNR at a stage never before observed, 
and to learn about the physics of shock acceleration in 
faster shocks than seen in any SNR (v s ~ 14, 000 km s" 1 ; 
Paper I). 

A simple model of a power-law electron distribution 
with an exponential cutoff at energy E max producing 
synchrotron radiation (XSPEC model sreut) has proved 
to be a useful tool in understanding X-ray synchrotron 
spectra in those doze n or so S NRs in which the phe- 
nomenon is observed (Reynolds 20081. The synchrotron 
spectrum cuts off more slowly than exponential, roughly 
as Sx exp[— (jz/j/ron) 1 / 2 ]. In addition to foreground 
absorption, this model requires three parameters: a 1 
GHz radio flux Sq, a mean radio-to- X-ray spectral index 
a (Sv oc v~ a ), and a rolloff frequency v m \\, the "criti- 
cal frequency" for electrons with energy E max , related to 
E ma , x by E 

max — 

39(/H/ ro u/l keV) 1 / 2 (B/10 [iG) 1 / 2 TeV. 
In applying this model in Paper I, we used Sg = 0.9 Jy, 
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Fig. 1. — Left: Chandra image of Gl.9+0.3, platelet smoothed 
(Willctt 2007). Colors are intensities only, between 1.5 and 6 keV. 
Image size 136" X 185". Right: 2008 radio image of Gl.9+0.3 
(Green et al. 2008): VLA at 4.9 GHz. Resolution 10" X 4". 

and obtained v ro n = 1.4 x 10 18 Hz, but with a very high 
absorbing column density, Njj — (5.5 + 0.3) x 10 22 cm~ 2 , 
implying significant scattering by dust along the line of 
sight. Such scattering removes photons from the source 
and distributes them in a faint halo out to arcminutc 
distances, but can also redistribute photons across the 
source, in both cases in an energy-dependent fashion. 

Below we present a reanalysis of the Chandra obser- 
vation using Markov chain Monte Carlo (MCMC) tech- 
niques, including the effects of dust scattering, with an 
emphasis on characterizing spatial variations. Such vari- 
ations hold important information on how the accelera- 
tion of electrons and /or magnetic-field amplification de- 
pends on different conditions such as shock speed and 
the obliquity angle O^n between the shock normal and 
upstream magnetic field, information crucial in under- 
standing shock acceleration in different astrophysical en- 
vironments. 

2. SPECTRAL ANALYSIS METHODS 

Radial averages of Chandra data show a dust-scattered 
halo out to about 3'. Thus, the spectral analysis of 
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Gl.9+0.3 must be conducted jointly with a halo anal- 
ysis. First, we refit the integrated spectrum without 
scatt ering, using a different a bundance set for absorp- 
tion (Grevesse & Sauval 19981, and a slightly different, 
improved background model. We then performed a joint 
analysis of the background-subtracted spectra from the 
source (innermost ellipse) and halo (between two outer 
ellipses) regions shown in Figure [2] 

We consider a simple uniform dust distribution along 
the line of sight to Gl.9+0.3, using tot al and differen - 
tial X-ray scattering cross sections from |Draine| ( |2003| . 
Small-angle scattering by dust out of the Gl.9+0.3 
source extraction region attenuates the X-ray spectrum 
by a factor of exp(-T sco ) + f src (l - exp(-r sco )), where 
r sca is the energy-dependent optical depth for scatter- 
ing, and f src is the energy-dependent fraction of pho- 
tons that have been scattered into the source extrac- 
tion region. Unless there are large quantities of dust 
in the immediate vicinity of the SNR, it is sufficient 
to consider only singly scattered photons in the second 
term containing f src . The fraction of si ngly scattered 



photo ns is equal to T sca exp(— r sca ) (e.g, Mathis & Lee 
19911, so scattering attenuates the source spectrum by 
(1 + fsrcTsca) exp(-T sca ). The fraction f src depends on 
the assumed dust distribution and the spatial structure 
of the X-ray emitting gas. For a point source, the angu- 
lar distribution of singly-sc attered photon s can be read- 
ily found from eq. (19) of Draine (2003). We approxi- 
mated the spatial structure of Gl.9+0.3 by the platelet- 
smoothed (Willett 2007) image shown in Figure [l] and 
convolved it with the model point source halo to arrive 
at f S rc- Similarly, the halo region flux is proportional to 
fhaio(l - exp(-r sca )), or /haioT sca exp(-T sca ) for singly 
scattered photons. We neglect multiply scattered pho- 
tons in modeling of the Gl.9+0.3 dust halo - this un- 
derpredicts the halo intensity at low energies where r sca 
exceeds unity, by at most 30% and probably less (Smith 
et al. 2006). We investigated the effect of the Chandra 
point-spread function (PSF), but found that only about 
1% of the SNR flux at high energies is scattered within 
the telescope into the halo extraction annulus. 
We performed joint modeling of the Gl.9+0.3 and halo 



spectra within the XSPEC software package (Arnaud 
1996 1, with dust scattering and PSF effects implemented 
as local multiplicative models. Model input parame- 
ters are energy-dependent functions f ps f (the energy- 
dependent fraction of the SNR flux scattered into the 
halo extraction region), f src , and fhaio, and a scalar pa- 
rameter Njf a that provides scaling for T sca in terms of 
the effective H column density for dust scattering. Un- 
der a single-scattering approximation, f src and fhaio de- 
pend only on the dust distribution along the line of sight, 
not on r sca , so they may be specified prior to modeling 
spectra. We use the sreut model described above for 
the intrinsic X-ray spectrum. While background sub- 
tracted spectra are shown for clarity in Figure [2] we do 
not subtract background prior to modeling so that we 
may use MCMC methods as implemented in XSPEC and 
PyMC 4 . We determined the normalization parameter Sg 
from our 1.5 GHz radio image, extrapolated to 1 GHz 

4 Patil, Huard, & Fonnesbeck (2009) PyMC: Markov 
chain Monte Carlo for P ython, version 2.0. Available at 



http:/ /pymc. googlecode.com 



using a spectral index of 0.62 ( Green et al. | 2008 1 , ob- 
taining a value of 1.17 Jy (under the possibly incorrect 
assumption that all the radio flux comes from the same 
population of electrons as produces the X-rays). The ab- 
sorbing and scattering column densities may differ if any 
of the assumptions underlying our uniform-dust model 
break down, e.g., an excess of dust local to the SNR will 
not affect absorption but will reduce large-angle scatter- 
ing, resulting in N H > Nff a . 

It is important to model the background correctly in 
view of the faintness of the dust scattered halo. We ex- 
tracted a background spectrum from a large (7.6' x 3.2' 
in size) region on the ACIS S3 chip beyond the dust- 
scattered halo. To separate sky and particle back- 
grounds, we also extracted spectra from the most recent 
stowed (particle only) ACIS S3 background file available 
online at the CXC 5 . This particle-only spectrum was fit 
with a combination of Gaussians and power laws, expo- 
nentially cut off on both ends. We then modeled our 
background spectrum as a combination of a sky back- 
ground, consisting of two absorbed power laws, and a 
particle-only background (allowing for change in inten- 
sity only). This fit, after appropriate scaling and again 
allowing for variations in the particle background inten- 
sity, was used in the modeling reported below. 

The Bayesian MCMC methods require specification of 
priors on unknown (fitted) parameters. Uniform, non- 
informative priors were used for Njj, Nff a , a, log^oii, 
r, and log Fx, with the following exceptions. In the 
spatially-resolved five-region analysis, we employed a 
normal prior for a with the measured radio index of 0.62 
and estimated standard (1 a) error of 0.02. A normal 
prior for Njj was assumed for protruding SE and NW 
extensions ("ears"), based on results of that fit. Since 
results are prior-dependent in this case, we also imposed 
an upper cutoff of v IO \i = 10 21 Hz on the corresponding 
prior, as theoretical considerations make higher rolloff 
frequencies unlikely. 

3. MODELING AND RESULTS 

We quantified the azimuthal brightness variations ob- 
vious in Figure [l] by circumscribing the shell with an 
elliptical annulus of width 19", and dividing it into 18 
equal segments. We found about a factor of 6 variation 
between the bright NW and SE limbs and the fainter 
regions between. 

Table [T] shows results for modeling the integrated spec- 
trum with the sreut model. Column 1 shows the im- 
proved no-scattering parameter values (means derived 
from MCMC chains) and their 5 and 95 percentile scores. 
Column 2 shows the corresponding values with our sim- 
ple scattering model. The chief difference is a lowering 
of v ro \\ by about 50%. The new value, 5.4 x 10 17 Hz 
(hv TO i\ = 2.2 keV) is still among the highest for a shell 
SNR. In the third and fourth columns we give results 
from a two-zone model, with the zones comprising most 
of the bright emission from the SE and NW limbs, in 
addition to a halo. Absorptions were tied together for 
the two regions, while the contribution from scattering 
from the "missing" faint limbs and center was accounted 
for by artificially increasing scattering by 20%, required 
to obtain consistency in the derived Nh values with the 

5 http:/ /cxc. harvard.edu/contrib/maxim/acisbg 
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value from the spatially-integrated model. The spectral 
indices a differ at the 99.98% level, and v ro \\ values at 
the 94% level. A similar simulation neglecting dust scat- 
tering gave qualitatively similar results but higher v 10 \\ 
values by over a factor of 2, indicating the importance of 
including dust scattering. These regions include about 
80% of the X-ray counts, but less than half the radio, 
suggesting that some of the radio flux may not be asso- 
ciated with the blast wave. 

Spectral results at higher angular resolution are desir- 
able. We subdivided Gl.9+0.3 into six regions as shown 
in Figure [2] but the X-ray data were not adequate to 
permit the use of the scattering model. Since scattering 
moves photons both out of the source and from brighter 
regions into fainter ones, the predominant effect is to 
smooth gradients, in both total brightness and spectrum, 
so these results should underestimate the true spatial 
variations. For a model-independent description of the 
observed data, we described the six regions simultane- 
ously with power laws, tying the absorption columns to- 
gether, but we also performed srcut modeling. Here the 
clearly discrepant radio properties of the "ears" made it 
necessary to impose a different prior on a (flat) than for 
the five regions within the source ellipse, so the "ears" 
modeling was performed separately. We also bounded 
the flat prior on log u IQ n by log f ro ii(max) = 21.0. Table [2] 
gives both power-law photon indices T (N(hv) cx [hv) 
ph cm -2 s _1 ) and srcut parameters. 

The systematic trends with brightness are apparent: 
the opposing bright limbs of the shell have harder spectra 
(either lower T or higher f ro ii) than the center, and the 
spectra soften also moving to the interior. The projecting 
fainter "ears" beyond the E and W peaks were combined 
for better statistics, and have the flattest spectra of all. 
The correlation of X-ray brightness with v^w strongly 
suggests a single cause. Our primary observational result 
is this systematic variation of spectral hardness, corre- 
lated with brightness. A complete dust-scattering model, 
including the effects of scattering back into source re- 
gions, applied to the spatially resolved data will require 
better X-ray data, but while the numerical values for v TO \\ 
will change somewhat, we do not expect these trends to 
be reversed. 

4. DISCUSSION 
We summarize our results as follows: 

1. Including dust scattering lowers v ro \\ by about 50% 
in the integrated spectrum, and by factors of more 
than 2 in the two-region fit, significantly lowering 
the inferred maximum electron energies. However, 
the froii values are still among the highest ever re- 
ported. The mean spectrum between radio and X- 
rays is significantly flatter in the SE than in radio 
alone for our two-zone model. 

2. Both X-ray brightness and spectral flatness or v m \\ 
show strong bilateral symmetry. As one moves 
away from the bright SE and NW limbs in azimuth 
or toward the center, the mean surface brightness 
drops by at least 6, while values of v m \\ drop by 
about an order of magnitude. 

3. The radio structure is dominated by a bright shell 
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Fig. 2. — Top: Raw Chandra image of Gl.9+0.3, showing source 
(cyan) and background (green) regions. The six regions used for 
spatially-resolved modeling are indicated (the "ears" are the sum 
of the two irregular regions, and the faint limbs and center com- 
prise the remainder of the source). Colors are intensities only (ct 
pixel - ), between 1.5 and 6 keV. The image has been smoothed 
with a Gaussian with <r = 3 pixels (1.5"). Bottom: Integrated spec- 
trum of Gl.9+0.3 and of its halo (top and bottom, respectively), 
extracted from regions shown above, after background subtrac- 
tion and adaptive binning (both done for display purposes only). 
Model: srcut modified by scattering (see text). For the halo spec- 
trum, contributions from scattering by dust and within the tele- 
scope are shown separately (upper and lower dotted curves, respec- 
tively). Two point sources within the dust-scattering halo region 
(between outermost two ellipses) were excluded. 



with somewhat smaller peak radius than the X-ray 
shell, and a single intensity maximum. 

4. The "ears" extending beyond the SE and NW rims 
in X-rays show the hardest spectra of all, and are 
not detected at radio wavelengths. 

The relatively poor correlation of radio and X-ray 
brightness around the shell, and the smaller radius of 
the peak of the radio shell, suggest that much of the ra- 
dio emission may not originate at the outer blast wave, 
but further inside, perhaps due to electron acceleration 
and/or turbulent magnetic-field amplification at the con- 
tact discontinuity between shocked ISM and shocked 
ejecta, as is probably the case in Cas A (e.g., Cowsik 
& Sarkar 1984). We focus on the X-ray structure and 
spectra. 
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In principle, the systematic azimuthal brightness vari- 
ations could be due to an inhomogeneous upstream 
medium that has the bilateral symmetry of Gl.9+0.3. 
This explanation requires higher densities in the SE and 
NW, where they would be expected to impede the ex- 
pansion. However, the remnant has a larger SE-NW 
than NE-SW diameter. It is equally unlikely to imagine 
systematically larger magnetic fields on opposite sides 
of the remnant. However, the observed simple bilateral 
symmetry of both brightness and spectral variations can 
arise naturally if electron injection, acceleration rate, or 
magnetic-field amplification depend somehow on shock 
obliquity #Bn- (The physical origin of such dependences 
is not clear.) If the magnetic field is roughly uniform 
and not too far from the plane of the sky, then as a 
spherical shock encounters it, the obliquity angle around 
the remnant periphery will vary with the observed bilat- 
eral symmetry. That explanation essentially requires the 
supernova to be Type la, as a core-collapse supernova 
is expected to expand into stellar-wind material with a 
frozen-in magnetic field in a tightly wound Parker spiral 
(roughly azimuthally oriented) with all obliquities there- 
fore near 90°. 

Synchrotron brightness from a power-law distribution 
of electrons N(E) = KE-( 2ct+ V varies as KB 1+a . The 
electron energy density can vary with obliquity if the in- 
jection or acceleration rate depends on #Bn, while the 
magnetic-field strength will vary due to simple flux- 
freezing and shock compression of the tangential com- 
ponent, and possibly also to turbulent a mplification due 
to co smic-ray-induced instabilities (e.g., |Bell & Lucek 



diffusion properties (E m3 ) (Reynolds 1998). These de- 
pend differently on physical parameters, including shock 



2001 1, a process which might also have obliquity depen- 



dence. However, in all known cases, X-ray synchrotron 
emission is produced by electrons on the loss-steepened 



tail of a p ower-law distribution (e.g., |Reynolds fc Keo- 
hane| |l999), introducing a third parameter. For the sim- 



plc srcut model, that parameter is the rolloff frequency 
f ro n- X-ray brightness will vary even for constant K and 
B if v TO \\ varies, and brighter regions will have harder 
spectra. 

For a uniform upstream magnetic field, if the post- 
shock field is simply compressed, a strong-shock com- 
pression ratio of 4 will produce post-shock magnetic 
fields four times larger where the shock is perpendicular 
($Bn^ 90°), and hence a larger synchrotron emissivity by 
4 1+Q = 9.5 compared to where #Bn~ 0° (using a = 0.62; 
Green et al. 2008). This is adequate to explain the re- 
quired factor of about 6, but would imply that the bright 
limbs are projections of an equatorial "belt" which might 
be detectable across the remnant center in a deeper ob- 
servation. 

Standard first-order Fermi shock acceleration in the 
test-particle limit predicts a spectral index a — 0.5 for 
a strong shock. Steeper radio indices are universal in 
young remnants, perhaps indicating cosmic-ray-modificd 
shocks (e.g., Reynolds & Ellison 1992). We find a steeper 
radio-to- X-ray value of a than that in the radio alone, 
consistent with the expectation of modified shocks. Bet- 
ter radio data will be required to confirm the radio index 
and search for possible spatial variations. 

In Fermi shock acceleration, three possible limitations 
on E max are synchrotron losses (E m i) 7 finite age (or size; 
E m 2), and escape due to an abrupt change in upstream 



obliquity, which may affect the acceleration rates ( Jokipi: 
1987). If we lump all such obliquity dependence into 
a factor ftbn(@Bn), the corresponding rolloff frequencies 
obey u IO \i cx £^ ax i? and we have 



V m \ OC ftbn u 8 v m2 cx f tbn B 2 t u s v m?> cx A m B x B 2 

(!) 

where B\ and B 2 are the upstream and downstream 
magnetic-field strengths, respectively, ug, = u s h/10 8 cm 
s _1 , and X m is a wavelength of MHD waves above which 
the diffusion coefficient is assumed to jump to much 
larger values, mimicking the absence of longer waves ca- 
pable of scattering more energetic particles. 

The absence of _B-dependence in v m \ means that if elec- 
tron energies are limited by radiative losses, the different 
spectral slopes cannot be produced purely by varying 
the magnetic field. The strong shock-speed dependence 
could produce a significant effect; but the observed order- 
of-magnitude variation in ^ ro n would require a greater 
variation in it s h than seems consistent with the remnant's 
circular outline to within about 10%. 

Since B 2 jB\ is obliquity-dependent, either age-limited 
(case 2) or escape-limited (case 3) acceleration can ex- 
plain some spectral variations with obliquity even if 
ftbn = 1- A non- linear variant of escape-limited accel- 
eration in which A m grows as increasingly higher-energy 
particles escape might be able to provide the observed 
range; otherwise one expects A m roughly constant, a 
property of the upstream medium, giving a maximum 
variation of v m3 of 4 or so. Age-limited acceleration 
could easily accomplish this without additional obliquity- 
dependence, since B\ will vary by a factor of at least 64. 
We conclude that a minimal explanation of both bright- 
ness and spectral variations is interaction of the blast 
wave with a roughly uniform magnetic field parallel to 
the bright rims (roughly SW - NE), with no additional 
obliquity-dependence of either electron injection or of ac- 
celeration rates. In that case, the "ears" must represent 
enhanced upstream diffusion where the shock is perpen- 
dicular, which seems unlikely since cross-field diffusion 
mean free paths are probably limited to the gyroradius. 
However, a preshock-diffusion explanation is consistent 
with the lack of radio emission from the "ears" (as lower- 
energy radio-emitting electrons would not diffuse as far 
ahead of the shock). In this case, an srcut model for the 
"ears" would not be appropriate, as already suggested by 
their discrepant value of a. 

An alternative explanation involves a magnetic field in 
which the shock is parallel at the "ears." This would 
almost certainly require enough turbulent amplification 
of magnetic field to swamp the frozen-in field increase 
where the shock is perpendicular. In addition, the mag- 
netic field would need to be quite close to the plane of 
the sky, so that the bright limbs, now "polar caps" seen 
edge-on, are at the edge of the remnant in projection. 
The same requirement occurs for this geometry in SN 
1006. Strong magnetic-field amplification probably re- 
quires efficient shock acceleration, so that a population of 
energetic ions is implied. However, for any magnetic-field 
geometry, effects due solely to magnetic-field dependence 
of acceleration time cannot produce the spectral varia- 
tions if electron acceleration is limited by synchrotron 
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TABLE 1 

Spectral Fits 



Region 


Whole 


Whole 


Two regions (scat) 




(no scat) 


(scat) 


SE NW 


N H (abs) (xlO 22 cm- 2 ) 


6.76 (6.37, 7.16) 


5.06 (4.48, 5.62) 


5.12 (4.56, 5.68) 


Nff (scat) 




3.48 (2.68, 4.35) 


3.60 (2.83, 4.41) 


S 9 (Jy) 


1.17 


1.17 


0.203 0.291 


a 


0.649 (0.625, 0.673) 


0.634 (0.614, 0.655) 


0.566 (0.541, 0.592) 0.612(0.588, 0.637) 


log i^roll (Hz) 


17.90 (17.58, 18.29) 


17.73 (17.48, 18.01) 


17.48 (17.23, 17.76) 17.72 (17.44, 18.06) 



Note. — Errors are 90% confidence limits throughout. The no-scattering fit updates the values from Paper I. 



TABLE 2 

Multi-Region Spectral Parameters 



Region 


Sg (mjy) 


a 


log I^roll (Hz) 


r 


Flux 


Outer SE 


67 


0.60 ± 0.02 


18.12 (17.80, 18.52) 


2.2 (2.1, 2.4) 


6.5 (6.2, 6.9) 


Inner SE 


190 


0.62 ± 0.02 


17.46 (17.24, 17.70) 


2.6 (2.4, 2.8) 


4.4 (4.2, 4.7) 


Faint limbs, center 


600 


0.63 ± 0.02 


17.16 (16.98, 17.35) 


2.8 (2.6, 3.0) 


5.8 (5.4, 6.2) 


Outer NW 


68 


0.63 ± 0.02 


18.55 (18.06, 19.21) 


2.1 (1.9, 2.3) 


4.8 (4.5, 5.1) 


Inner NW 


220 


0.64 ±0.03 


17.69 (17.43, 18.00) 


2.4 (2.2, 2.6) 


4.7 (4.4, 5.0) 


"Ears" 


< 2 


< 0.51 ± 0.04 


< 18.67 (17.65, 20.37) 


2.0 (1.7, 2.2) 


3.4 (3.1, 3.7) 



Note. — Errors are 90% confidence limits. Fitted absorption N H = 7.3 (6.8, 7.7) X 10 22 cm 2 for power-law models and 6.9 (6.6, 7.1) X 
10 22 cm -2 for sreut models. Flux is in units of 10 — 13 ergs cm -2 s~ 1 between 1 and 8 keV. "Ears" sreut modeling was performed 



independently: we assumed a flat prior on v TO \\ to a maximum log i/ T , 

losses. 

Currently scheduled observations may allow consider- 
able clarification. Our 250-ks Chandra follow-up, and 
the higher- resolution VLA study we have undertaken, 
should allow better understanding of the relation of radio 
to X-ray electron populations, better characterization of 
the scattering, and better deduction of spatially resolved 



: 21.0, and a flat prior on a. 

spectral properties. Gl.9+0.3 will provide valuable con- 
straints on the obliquity-dependence of shock accelera- 
tion and magnetic-field amplification. 
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